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We report on our theoretical study of the magnetoresistance in spin polarized transport through 
a finite carbon nanotube (CNT). Varying the Fermi energy of a CNT and the relative strength of 
couplings to two ferromagnetic (FM) electrodes, we studied the conductance as well as the magne- 
toresistance (MR). Due to resonant transport through discrete energy levels in a finite CNT, the 
conductance and MR are oscillating as a function of the CNT Fermi energy. The MR is peaked at 
the conductance valleys and dipped close to the conductance peaks. When couplings to two FM 
electrodes are asymmetric, the MR dips become negative under a rather strong asymmetry. When 
couplings are more or less symmetric, the MR dips remain positive except for a very strong coupling 
case. Under strong coupling case, the line broadening is significant and transport channels through 
neighboring energy levels in a CNT interfere with each other, leading to the negative MR. 

PACS numbers: 72.25.-b, 73.40.Gk 



I. INTRODUCTION 

Spin-polarized transport^ has attracted lots of atten- 
tion because of its potential applications to spintronic de- 
vices. Spintronics attempts to use the electron spin in or- 
der to control the electric current in the devices. The spin 
valve, consisting of a non-magnetic (NM) material sand- 
wiched in between two ferromagnetic (FM) electrodes, is 
a typical two-terminal spintronic device. The current in 
the spin valve is modulated by the relative orientation of 
magnetizations in the two FM electrodes. Usually more 
current flows when the two magnetizations are parallel 
than when they are antiparallcl. The difference in re- 
sistance between two configurations of magnetization is 
called the magnetoresistance (MR). 

In addition to the relative orientation of magnetiza- 
tions, there are many other methods to control the cur- 
rent in the spin valve systems. In magnetic tunnel junc- 
tions (MTJs), the oxide barriers can change the tun- 
neling magnetoresistance (TMR) in a significant way. 
With MgO barrier, TMR values reaching several hun- 
dred percents were predicted theoretically 2 - and reported 
experimentally. 3 The MgO layer acts as the highly spin 
selective filter. The sign of TMR can be changed with 
different oxide barriers. For example, tunneling through 
Al oxide barriers is dominated by s electrons leading to 
normal positive TMR. On the other hand, the d electrons 
can tunnel SrTi03 barriers 4 more easily with the negative 
or inverse TMR. The oxidation process of the insulating 
barriers also affects the TMR values. Strongly temper- 
ature dependent suppression of TMR was reported^ for 
MTJs with heavily oxidized Al oxide barriers. Inclusion 
of the dusting layers^ in the MTJs changes the TMR 
values, too. 

The current in the spin valve can also be modulated 
by the third electrode which is coupled capacitively to 
the non-magnetic material. Though the three-terminal 



spintronic device or the spin field-effect transistor (FET) 
was suggested theoretically^ it was not realized experi- 
mentally as yet. To implement the gate electrode into 
the spin valve tends to lengthen naturally the nonmag- 
netic part. The success of the spin FET seems to depend 
on the efficient spin injection from the FM electrode into 
the NM part and the preservation of the spin coherence 
in the NM part. Due to the low atomic number of car- 
bon and accordingly the weak spin-orbit interaction in 
carbon systems, the CNT and graphene are considered 
as ideal candidates for the NM part in the spin valves. 
The CNTs and graphenes are featured with a very long 
spin-diffusion length (up to a few fim) and spin-flip time 
(up to a few tens of ns) at room temperature. 

There have been several experimental attempts 
to measure the MR in the spin valves with the 
CNTs^&^dmdldlds the graphenes^^^ and the 
fullerenes (C60)^i2i For the spin valves with the CNT, 
the positive MR&i^ as well as negative MRi 3 - were re- 
ported. The inverse MR (up to 6 %) was observed 13 
in the Co/CNT/Co and Co/CNT/Ni spin valves with a 
highly transmissive contact. On the other hand, the pos- 
itive 61 % MR is observed^ 2 - in the CNT spin valves with 




FIG. 1: Schematic display of a carbon nanotube (CNT) sys- 
tem. The finite CNT consists of TV layers and is end-contacted 
to the left and right ferromagnetic electrodes. 
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the half- metal electrode (tunneling contact). More in- 
terestingly, some experimental group a 14 i 15 succeeded in 
modulating the MR in the CNT spin valves (NiPd elec- 
trodes) with the gate electrode. The observed MR and 
conductance are correlated and oscillating as a function 
of the gate voltage. The MR is positive and peaked in 
the conductance valleys, while it is suppressed and be- 
comes even negative near the conductance peaks. The 
control of the MR and conductance with the gate volt- 
age suggests that the resonant tunneling through quan- 
tized energy levels in a finite CNT is responsible for the 
spin-polarized transport. 

Some experimental features can be understood in 
terms of the spin-polarized transport through a single 
resonant levels Suppose that one energy level with en- 
ergy eo is coupled to two source and drain (left and right) 
ferromagnetic electrodes. The transmission probability 
for an electron with spin direction a — ± (spin- up/spin- 
down) from one electrode to the other is given by the 
expression at the Fermi level in the linear regime 



T o (0) 



el + (r La + r Ra y 
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eo can be adjusted by the gate voltage. T pa is the spin- 
dependent linewidth of a resonant level coming from the 
coupling to the p — L,R (left and right) electrode. The 
magnetoresistance is dependent sensitively on the posi- 
tion of eo relative to the Fermi level in the electrodes. 
Close to or on resonance (eo = 0), the transmission prob- 
ability is simplified as 
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For highly asymmetric coupling to two electrodes, say, 
Tl r^, the transmission probability is simplified as 
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and the MR ratio can be obtained as 

Gp - Gap 20l(3r 



MR = 
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The inverse or negative MR is obtained for a highly asym- 
metric coupling case. Here Gp and Gap are the conduc- 
tance in the spin valve for the parallel and antiparallcl 
configurations of magnetizations, respectively. The spin 
polarization (3 P in the linewidth is defined by the relation 
T p ± = T p (l ± f3 p ), where p — L,R denotes the left and 
right electrode. 

Off resonance or when the Fermi level is far away from 
the resonant level, the transmission probability can be 
simplified as 



^La^Ra 



(5) 



This approximation gets better with the smaller 
linewidth. The MR ratio can be readily computed from 
this approximation. 
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This expression of MR is similar to that for the Julliere 
model of the magnetic tunnel junction (MTJ). From this 
simple resonant level model, we can deduce that the MR 
is bounded roughly by two values of Eqs. f4| and ([6|). The 
MR is dipped at the value of Eq. ((4J near the conductance 
peaks and is peaked at the value of Eq. © at the con- 
ductance valleys. Though some features of experimental 
MR can be understood based on a single resonant level 
model, we shall show in this work that the interference 
between multi quantized energy levels in a finite CNT 
plays an important role in the spin-polarized transport. 

In this paper we study theoretically the magnetore- 
sistance in the phase-coherent spin-polarized transport 
through a finite carbon nanotube (CNT). Our model sys- 
tem, schematically shown in Fig.[TJ is a typical spin valve 
with a carbon nanotube sandwiched in between two (left 
and right) FM electrodes. In addition, the carbon nan- 
otube is coupled capacitively to a gate electrode such 
that the energy levels or the Fermi energy in a CNT 
can be shifted up and down by the gate voltage. In our 
model study, the control parameters are the Fermi en- 
ergy level in a CNT (or the gate voltage) as well as the 
relative strength of coupling constants between a CNT 
and two FM electrodes. The linear conductance as well 
as MR are oscillating as a function of the CNT Fermi en- 
ergy. The MR is featured with a peak in the conductance 
valleys and a dip structure near the conductance peaks. 
The shape of MR as a function of the gate voltage de- 
pends on the relative magnitude of couplings to the left 
and right FM electrodes. (1) With asymmetric couplings 
to the two FM electrodes, the MR dips become nega- 
tive under a high asymmetric aspect ratio of couplings. 
When the coupling strength is increased, the (negative) 
MR dip structure is broadened. (2) For symmetric cou- 
plings, qualitatively different behavior is observed in MR 
depending on the coupling strength. In the case of weak 
couplings, the MR as a function of the gate voltage is 
oscillating without any sign change. The simple peak 
appears in the conductance valleys, but the positive dip 
near the conductance peaks has an additional local peak, 
that is, a dip-peak-dip structure. When the coupling 
strength is increased, discrete energy levels in a CNT 
are broadened and overlap each other. Due to the inter- 
ference between neighboring energy levels, the dip now 
becomes negative and the MR shape is highly asymmet- 
ric with respect to the peak position. One remarkable 
point is that the inverse MR in resonant transport can 
be observed for the case of symmetric strong couplings or 
for the highly transmissive contact between CNT and the 
FM electrode. Our study suggests another way to change 
the MR values by controlling the coupling between FM 
electrodes and a finite CNT. Our study can be equally ap- 
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plied to other nanostructures where discrete energy levels 
are formed due to their finite size. The preliminary re- 
sults of our work is already reported elsewhere^ 

The rest of this paper is organized as follows. In Sec. 
II, the model Hamiltonian is introduced for the FM- 
CNT-FM system and the spin-polarized current is for- 
mulated in a Landauer-Biittiker form. The results of our 
work are presented in Sec. Ill and a conclusion is in- 
cluded in Sec. IV. In Appendix [M we present the tight- 
binding Hamiltonian approach to a CNT and the phase 
information about the tt and tt* states. In Appendices 
|B1 and [Cj we elaborate on the coupling matrix between 
the CNT and the ferromagnetic electrodes and discuss 
its symmetry. 

II. FORMALISM 

To study the phase-coherent spin-polarized transport 
through a carbon nanotube, we consider the model sys- 
tem which is schematically displayed in Fig. Q] The fi- 
nite armchair-type (n, n) carbon nanotube (CNT) is end- 
contacted to the two ferromagnetic electrodes. The band 
structure of the metallic CNTs close to the Fermi level^i 
is known to be described accurately by the single tt elec- 
tron tight-binding Hamiltonian 

H cnt = e 9 a L a *a - 1 X! ( a L a 3c* + h - c ) 

ia <zj> a 

= (7) 

Here i, j runs through the atomic carbon sites in CNT 
and < i,j > denotes the nearest neighbor pairs, a — ± 
represents two spin directions, up (+) and down (— ), and 
the hopping integral t = 2.66 eV^i The on-site energy e g 
at each carbon site is proportional to the gate bias voltage 
and is chosen to be zero in the absence of the gate bias. 
Strictly speaking, the on-site energy at each carbon site 
depends on the separation of each site from the substrate 
(which is capacitively coupled to the gate electrode). In 
this work, we ignore such dependence and assume the 
on-site energy to shift uniformly by the gate voltage. In 
the second line, H cnt is written in a matrix form, where 
H C nt is the square matrix and reflects the structure of the 
tight-binding Hamiltonian in the (n, n) CNT. ty a is the 
electron annihilation operator represented by the column 
vector and has as many components as the number of 
atomic sites in the finite CNT. 

The (n, n) CNT has the discrete rotational axial sym- 
metry C nv about its axis. For details, look at Appen- 
dices [3] and The energy bands in a CNT can be clas- 
sified using the group representation. Since the energy 
dispersion is a function of wave vector along the nan- 
otube axis, every energy band can be uniquely specified 
by the quantum number belonging to the irreducible rep- 
resentations. There are in total An bands for the (n, n) 
carbon nanotubes out of which four bands are nondegen- 
erate and other 2n — 2 bands are doubly degenerate. The 



nondegenerate four bands are named as tt and tt* bands 
(two tt and two tt* bands). Out of these four bands, only 
two (tt and tt* bands) cross the Fermi level (Ep). Only 
these tt and tt* bands, crossing the Fermi level, are rel- 
evant to our spin polarized transport study. All other 
bands are gapped close to the Fermi energy. 

In finite (n, n) CNTs with N layers, the discrete energy 
levels close to Ep of tt (bonding) and tt* (antibonding) 
characters will provide the transport channels for elec- 
trons incident from the ferromagnetic electrode. Out- 
side this energy window, other higher-lying bands start 
to contribute to the electronic transport. The tt states 
are featured with no change in the relative phase from 
one carbon site to the other along the circumference in 
one layer. On the other hand, the tt* states have the 
alternating relative phase, +1 and —1. For details, look 
at Appendix lAl 

The left (p = L) and right (p = R) FM metals are 
described by the two conduction bands of majority and 
minority spins. 

H <* = E J2 e pk a c U a c P k a - (8) 

P=L,R k a 

Here c' - and c r are the creation and annihilation op- 

pka pka r 

erators, respectively, for electrons of wave vector k in 
the electrode p = L,R, with the spin direction a = ±. 
a = +(— ) means the spin is aligned parallel (antiparallel) 
to the direction of magnetization. e p £ a is the energy dis- 
persion relation for electrons in the ferromagnetic metals. 

The coupling between a CNT and the electrodes is nei- 
ther well controlled experimentally nor well understood 
theoretically yet. In one theoretical paper— only the tt* 
states are coupled to the electrode of a jellium model 
and the tt states are effectively decoupled. In other the- 
oretical work)2& the tt (tt*) states are strongly (weakly) 
coupled to the Al electrodes for the end-contact geome- 
try. In this paper we are going to adopt the symmetry- 
adapted coupling based on the group theory and confine 
our interest to one specific end-contact geometry. For 
the end-contact geometry, only the carbon atoms at the 
left and right edge layers are assumed to be coupled to 
the FM electrodes. This coupling Hamiltonian can be 
written as2£ 

pka 3 

V is the coupling Hamiltonian between the CNT and 
FM electrodes and assumed to be spin-conserving. The 
coupling can be dependent on the electron spin direction. 

In order to find the electric current we need to rewrite 
the model Hamiltonian in the symmetry- adapted basis. 
From here on, without loss of any generality we focus 
our discussion on the (5, 5) CNT. We can define the an- 
gular momentum m about the nanotube axis for elec- 
tronic states in a CNT as well as in the FM electrodes. 
However the angular momentum states e lm ^ (m is an in- 
teger) are not the basis functions of the symmetry group 
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of the (5, 5) CNT. As discussed in Appendix [Cl the tt 
(tt*) states belong to the Ai (A 2 ) irreducible representa- 
tion (irred. rep.), respectively. The basis functions of A\ 
irred. rep. are cos(5m0) with nonnegative integer m and 
those of Ai irred. rep. are sin(5m</>) with positive inte- 
ger m. The spin polarized electrons should have the wave 
functions with angular dependence of cos 5m</> (sin 5m<j)) 
in order to tunnel resonantly through the tt (tt*) states 
of the (5, 5) CNT, respectively. 

We project the electron states of the ferromagnetic 
electrodes based on the irred. reps, of the (5, 5) CNT 
and keep only the states belonging to the A\ and A 2 ir- 
red. reps. Other projected states belonging to E\ and 
Ei representations are not coupled to the tt and tt* states 
of the (5, 5) CNT. 

Hi = J2J2cK(pka\V\ja)a ja +H.c. 



Here (f>j represents the angular position of the j-th car- 
bon atom along the circumference of the edge layer. The 
angle is measured from any line which intersects the tube 
axis perpendicularly and bisects two neighboring carbon 
atoms at the edge layer. For example, 4>j takes the fol- 
lowing set of values for (5, 5) CNT, alternating between 
the odd-th and even-th layers. 



5 \ 6' 6' 6 6' 6 

2tt f 1 1 1 1 1 . 
— \ ±-,l± -,2± -,3± -,4± - \ . 12) 
5 13' 3' 3' 3' 3 ' v ' 



We need to consider only those angular momentum states 
m which are multiples of 5 for the (5, 5) CNT system. 



k P3 a 



kzk± m= — oo pja 

+H.c. (10) 

Here we note that k — k z z + k± and k± — \k±\. The az- 
imuthal angle dependence can be written down explicitly 
as 

(pk z k±ma\V\ja) = V pma (k z , k x )e im ^ . (11) 



H * = X X Y, c kk ±5ma V p5ma (k z ,k^)e^a ja 

k^k± m= — oo pja 

+H.c. (13) 



Rearranging terms in the basis functions of the A\ and Ai 
irreducible representations, the relevant coupling Hamil- 
tonian can be written down as H\ = + H\-k* ■ 



H\ v 

Hlyr* 



oo 



k z k± m—0 pja 



XXX [ C ^*fc,fc imQ V*ma( fc ^^) sin ( 5m ^) a ja + H -C. 
k z kj_ m—1 pja 



(14a) 
(14b) 



In our model study, two projected tt and tt* bands con- 
tribute to the electronic transport independently. Each 
band has mutually independent infinite number of trans- 
port channels (indexed by m). For angular position <pj 
given by Eq. (|12[) . cos(5m<j)j) is independent of j while 
sin(5?7j0j) oc (— iy . That is to say, the projected tt and 
tt* conduction bands correctly pick up the tt and tt* states 
in the CNT, respectively. 

According to the above analysis based upon the sym- 
metry, the tt and tt* states are characterized by the fol- 
lowing properties: (i) All carbon atoms at the edge layer 
are coupled to the electrode with the same magnitude of 
coupling constants; (ii) The coupling constants are uni- 
form in phase for the tt states, but alternating in sign 
from one atom to the other for the tt* states. These 
properties are true for all projected bands in the elec- 
trode. We can introduce the following effective model 



Hamiltonian 

p—L,R eja 

+ X T, c U*eaV P ^a(e)(-l) j a ja +H.c{15) 

p=L,R eja 

Here e is the energy of electrons in the ferromagnetic 
electrode. This effective Hamiltonian is totally equivalent 
to Eq. fT4")) , as far as the transport is concerned. Two 
different tt and tt* conduction bands are coupled to the 
CNT with the effective coupling constants defined by the 
relation 

Np Va \V p7ra \ 2 = ^ Npirmg | Vpyma | 2 COS 2 {mTT/3) ,( 16) 
m 

N p7r *a\V P 7r*a\ 2 = ^ N p7r *ma \ Vpir* ma | 2 sill 2 (m7r/3Ql7) 

m 
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Here N and V denote the density of states and coupling 
strength, respectively. The coupling in the n band is in- 
dependent of j or is uniform in phase, while the coupling 
in the ir* band is alternating in its phase. Due to this 
property in the coupling, the tt and tt* bands are correctly 
coupled to the tt and ir* states in the CNT, respectively. 

In this paper we are going to consider only the collinear 
magnetizations or parallel and antiparallel alignment of 
magnetizations in the two FM electrodes. Using the 
nonequilibrium Green's function method ) 27 ' 28 ! 29 we can 
readily derive the spin-polarized current flowing from the 
left electrode to the right one and it can be written as 

30 .31.32 

I La = \ J deT a (e) [f R ( e )-f L (e)}, (18a) 
T a {e) = 4TrD r a r La D«T Ra . (18b) 

Here /l/r is the Fermi-Dirac distribution function in 

the left/right electrode. D r J a is the retarded/ advanced 
Green's functions of the CNT whose self-energy is deter- 
mined by excursion of electrons in the CNT into the two 
ferromagnetic electrodes and is given by the expression 

K /a (t) = ^V po (t)(^'(fc l£ )^(k). (19) 

pk 

Here V p is the coupling represented in a column vec- 
tor and has 2n components. For example, VJ a = 
V p7Ta (l 7 1, 1, 1, • ■ ■ , 1) for the tt states and V^ a — 

VjrK*a{— 1) 1) — 1) 1) • • • ,1) for the states. G r Ja is the 
retarded/ advanced Green's function of electrons with 
spin direction a = ± in the electrode p = L, R. Hence 
the Green's function of the CNT can be written as 

D r a (e) = [el - Hot- E;(e)]- X . (20) 

The hybridization matrices T pa = — ImSp Q are the imag- 
inary part of the self-energy 

r pQ (e) - 7rJ2V pa (k)VX(k) 6(e-e pka ). (21) 

k 

Assuming the constant coupling constant we have 

T pa = TrN pa {e)V pa Vj a . (22) 

In components, T pa ^ = tt N pa V pal V* aj . Here N pa (e) 
is the density of states for the p = L,R FM electrode 
for spin direction a = ±. For the tt states, r p7rQ .y = 

^■^piza. I ^p7ra I and ^ P 7T*a,ij ( 1) ^^"^Jj7r*a|^ / p7r*a| ft->r 

the tt* states. The expression of spin-polarized current in 
this paragraph applies to the projected tt and tt* bands 
separately. The total current is the sum of two contribu- 
tions. The hybridization or lincwidth T p ± is parameter- 
ized as 

r P ± = ttN p± \v p± \ 2 = r p (i±p p ) (23) 

Here fi p is the effective spin polarization of the FM elec- 
trode p = L,R and its value is chosen as (3l = = 0.2 
in our numerical works. 



III. RESULTS AND DISCUSSION 

We computed the magnetoresistance (MR) as well as 
the linear response conductance as a function of the on- 
site energy at the carbon atomic sites while varying the 
diameter (n) and the length (N: the number of layers) of 
the CNTs. The effective length of a CNT with N layers 
is L = Na/2 with a lattice constant a = 2.46 A. Though 
the details of the results depend on both n and N, we 
can extract out the generic features in the conductance 
and MR. For the presentation, we have chosen the (5, 5) 
carbon nanotubc (n = 5) with N = 999, 1000, 1001 lay- 
ers. Using the property of repeating layer structure, the 
desired Green's function is computed in a recursive way 
and the relevant matrix dimension is reduced to 2n x 2n. 

The generic features in MR and conductance do not 
show any obvious even-odd parity effect of N close to 
the Fermi level (zero gate voltage), but instead depend 
on N in modulo 3 as displayed in Figs. [2j [3l and |U 
For example, the results for N = 999 and N = 1002 
are similar in their structure. The discrete energy level 
spacing decreases with increasing number of layers. For 
the end-contact geometry we used the symmetry-adapted 
coupling between the electrodes and carbon atoms at the 
edge layers. The MR is oscillating as a function of the 
gate voltage and has a dip in MR near the conductance 
peak. Furthermore the dip can become negative depend- 
ing on the asymmetry of the couplings to two ferromag- 
netic electrodes and on the strength of the couplings. 

According to the recent ab initio calculations , 26 ' 33 the 
coupling strength between carbon atoms in the CNT and 
the electrodes depends on the atomic elements in the elec- 
trode. For example, Au and Al are relatively weakly cou- 
pled to the carbon atoms, while Ti electrodes are strongly 
coupled. The Pd electrode is known to be most strongly 
coupled to the carbon nanotubesi 34 ' 35 Furthermore, the 
coupling strength varies from sample to sample and is de- 
pendent on the sample fabrication process. For our model 
study, we treat the strength of the coupling constants as 
adjustable model parameters in order to investigate its 
effect on the spin-polarized transport through the CNTs. 

To obtain some insight about the structure of coupling 
matrix elements between the ferromagnetic electrode and 
the finite (n, n) CNTs, we consider the simple jellium 
model— i^ 6 . for the ferromagnetic electrodes. As noted in 
the above, only the discrete energy levels of tt and 7r* 
states contribute to the electron transport close to the 
Fermi level. Accordingly, the conducting electronic states 
in the FM electrodes should have the same symmetry 
as the 7r and tt* states. As a simple estimation we may 
consider the plane wave on the ferromagnetic surface and 
expand it in terms of angular momentum states 

oo 

^k-x = e rk zZ i m J m {ki_ry m ^. (24) 

m— — oo 

Here x is the position vector on the ferromagnetic inter- 
face from the center of a carbon nanotube. k = k±_ + k z z 
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is the Fermi wave vector, and J m is the Bessel function. 
For (n, n) CNT, the tt and tt* states are characterized 
with the azimuthal angular momentum quantum num- 
ber being an integer multiple of n. Since the carbon 
atoms at the edge layer is more or less localized along 
the circumference of radius r n [r n is the radius of the 
(n, n) CNT], the coupling strength of the tt and tt* states 
can be roughly proportional to J n m(k±r n ). The m = 
term contributes only to the coupling of the tt state. 
Since the Bessel function is oscillating with its argument 
and its amplitude is proportional to l/>/k±r n for large 
k±r n (^> 1), the coupling constants may well depend on 
the Fermi wave number (kp) of the ferromagnetic elec- 
trodes and the radius [n and r n ) of a finite CNT. 

We cannot determine from our phenomenological ap- 
proach the coupling matrix elements, but instead use the 
group theory to find the symmetry-adapted coupling ma- 
trix elements. Considering the Jellium model for the 
ferromagnetic electrodes, we assume that the electronic 
density is relatively uniform on the FM interface. This 
will be true for the conducting s and p electrons^! On 
the other hand, the conducting d electrons have more 
or less a localized character such that the spatial varia- 
tion of the electron density cannot be considered uniform. 
Which electrons of the FM electrode, s and p or d, play a 
dominant role in the spin-polarized transport depend on 
the materials of the spintronic devices. As an example, 
the s and p electrons are main carriers in the magnetic 
tunnel junctions with the Al oxide barrier but the d elec- 
trons are responsible for the spin-polarized transport in 
the MTJ with the SrTiC>3 barrier. Since the coupling 
between the FM electrodes and a CNT is not known at 
present, we consider the case that the s and p electrons 
are more strongly coupled to the CNT. In this case each 
C atom at the edge layer is coupled to the FM electrode 
with the same strength. Even in this case, however, the 
relative strength of coupling constants for the tt and tt* 
states will depend on the radius of the CNT as well as 
the Fermi wave number of the FM electrode, as can be 
expected from the simple estimation from the Eq. (f2"l| . 
The nature of the transporting carriers, tt or tt* states, 
may be probed by changing the radius or n of the CNT. 

Based upon the axial symmetry of the (5, 5) CNT, the 
coupling matrix elements can be expanded in terms of 
the angular momentum states as2£ 
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FIG. 2: Oscillating MR and conductance with asymmetric 
couplings to two ferromagnetic electrodes. The number of 
layers N = 1001. The linewidth parameters Tl.r are mea- 
sured in unit of t. (a) Weak coupling case, (b) Strong cou- 
pling case. MR is peaked in the conductance valleys and is 
dipped near the conductance peaks. For sufficiently asym- 
metric couplings between two electrodes, the MR is negative 
in the dipped regions. 



< ka\V\ja > = V 0a + V ma e" n ^. (25) 

Here <pj represents the azimuthal angular position of car- 
bon atoms at the edge layer around the CNT circumfer- 
ence. The electrode index p = L, R is suppressed for the 
notational convenience. As shown in Appendix [Cl the tt 
and tt* states in a (5, 5) CNT are characterized with the 
angular momentum quantum number m = 0, 5, 10, 
The most general form of the coupling matrix elements 



for the tt and tt* states can be written as 

< ka\V„\ja > = Vnoa + Ksa cos 5<f>j -\ ,(26a) 

< ka\Vir* \ja > — Vir*5a sin5^j + • • • . (26b) 

As already noted in Sec.|TTl the coupling constant for the 
tt state does not depend on the index j but the coupling 
constant for the tt* state is proportional to (— l)- 7 . In our 
work we take into consideration the azimuthal angular 
dependence based on the effective Hamiltonian, Eq. (| 1 5|) . 
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A. Transport through tt channel or it* channel 



The Hamiltonian of CNT, Eq. ([7]), is invariant under 
the particle- hole transformation ai a — > (— l) l a\ a when 
e g = 0. This particle-hole symmetry means that the dis- 
crete energy levels for a finite CNT appear symmetrically 
with respect to zero energy. Furthermore, the tt* levels 
are exactly the particle-hole image of the tt levels. When 
the number of layers satisfies the relation N — 3fc — 1 
(k is a positive integer), both levels of tt and tt* states 
lie very close to zero energy in a particle-hole symmetric 
way. Due to the presence of these two levels, the lin- 
ear response conductance, when no gate bias is applied, 
is systematically higher when N = 3k — 1 than when 
N = 3k, 3k + 1. These properties are well explained in 
the conductance panels of Figs. [21 [3l and 01 The tt chan- 
nel is chosen for our presentation in this subsection. The 
tt* channel gives qualitatively the same results as the tt 
channel. 

In Fig. [21 the MR as well as the linear conductance are 
presented for the asymmetric couplings (with the asym- 
metric aspect ratio Tl/Tr = 5) between two FM elec- 
trodes and the CNT. The MR as well as the linear con- 
ductance show the oscillating behavior as a function of 
the Fermi energy of CNT which is proportional to the 
gate voltage. The MR is positive in the conductance val- 
leys, but is suppressed and becomes even negative near 
the conductance peaks. The negative MR is realized for 
relatively high asymmetric aspect ratio (approximately 
larger than 4). In Fig. [2fb), the coupling strengths are 
increased compared to Fig.[5Ja) with the same asymmet- 
ric aspect ratio. The conductance peaks overlap each 
other and the MR is oscillating between positive and 
negative, with an increased range of negative MR. The 
negative MR in some spin valve systems was already ex- 
plained in terms of the spin-polarized resonant tunnel- 
ing with the asymmetric couplings. 22 With the higher 
asymmetric aspect ratio Tl/Tr, the inverse MR is fur- 
ther increased— but is limited by the lower bound value 
-2(3 L p R /(l + (3 L (3 R ) [Eq. ©]. 

In the conductance valleys, the MR is maximized 
and its value is limited by the spin polarization of the 
linewidth or the hybridization constant. The upper 
bound value is given by Eq. The maximum MR 

is well obeyed even for multiple discrete energy levels, 
as far as the energy level spacing is much larger than 
the linewidth. In order to increase the MR ratio, the 
larger spin polarization (3 in the linewidth is essential. 
The value of (3 can be increased by either using the half- 
metallic electrode (100 % spin polarization) or using the 
highly spin-selective coupling between the CNT and the 
FM electrodes. Recently, the higher MR ratio in the 
CNT spin valve was reported^ using the half-metallic 
electrodes. 

In Fig. [3l the MR and the conductance are displayed 
for the symmetric couplings between the CNT and two 
ferromagnetic electrodes. The MR oscillates with e g or 
the gate voltage and is dipped with positive values for 
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FIG. 3: Oscillating MR and conductance with symmetric cou- 
plings to two ferromagnetic electrodes. N = 1000. (a) Weak 
coupling case: MR is peaked in the conductance valleys and 
is dipped near the conductance peaks with the dip-peak-dip 
structure. The MR always remains to be positive in the weak 
coupling case, (b) Strong coupling case: When the coupling 
becomes strong enough for neighboring conductance peaks to 
overlap, the interference leads to highly asymmetric shape of 
MR and negative values of MR. 



the weak coupling case near the conductance peaks. The 
suppressed MR has a dip-peak-dip structure, which is dif- 
ferent from a simple dip for the asymmetric case. When 
the coupling strength is increased, the MR dip becomes 
negative and the MR shape is highly asymmetric. This 
inverse MR in the strong symmetric coupling case is rem- 
iniscent of the same inverse MR in the spin polarized 
transport through a quantum point contact i 39 i In this 
case, the inverse MR is also realized when the transmis- 
sion probability for both spin directions is close to a unity. 
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B. Transport through n and n* channels 
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FIG. 4: Effect of tt and tt* channels on MR. N = 999. Con- 
ductance and MR are presented when both tt and tt* channels 
are responsible for transport. We have chosen the linewidth 
parameters as I\ = T^* . 



In spin polarized transport through a single resonant 
energy level or widely spaced multi resonant energy lev- 
els, the inverse MR is possible only for the asymmetric 
couplings. 22 The negative MR in the symmetric strong 
coupling case is a direct consequence of the interference 
between neighboring energy levels or conductance peaks. 
This inverse MR is also related to the one-dimensional 
structure of a carbon nanotube. This point will be fur- 
ther discussed in Sec. IIII CI 



Atoms in the FM electrodes, coupled to the CNT, 
will rearrange their atomic positions to accommodate the 
coupling more efficiently. Atoms may try to conform to 
the local symmetry of the CNT. Though the m = (s- 
wave) coupling strength (V) is dominant, other m =/= 
contribution or the azimuthal angle dependence of V may 
be present. We are going to address the effect of this issue 
on the spin polarized transport through the CNT. 

The azimuthal angle dependence of the coupling ma- 
trix means the existence of nonvanishing components of 
higher angular momentum states. According to Eq. (|26[) . 
we have nonvanishing V^5 Q , V^»5 Q and higher coupling 
components. As proved in Sec. [Ill this situation can 
be treated by considering the effective coupling Hamil- 
tonian, Eq. ([15]). 



P7T ( 1 i Pp ) , 



(27a) 
(27b) 



We assume the same spin polarization for both tt and 
tt* projected bands. For our numerical presentation we 
have chosen the hybridization or linewidth parameter as 
Tp7r = Tpff* or the equal coupling strength for the 7r and 
tt* states. For realistic contacts, two coupling constants 
may well be different. 

In Fig. 2] (asymmetric aspect ratio T^/Tr — 10), 
the MR and conductance are presented when both tt 
and tt* states in the CNT are responsible for the trans- 
port. As already mentioned in Sec. IIII Al the tt* lev- 
els are the particle-hole images of the tt levels for a fi- 
nite (n, n) CNT with our tight-binding Hamiltonian ap- 
proach, Eq. 0. The tt (tt*) states are located in the 
energy range —it < E < t {—t < E < 3i), respectively. 
As shown in Fig. [2J the conductance peaks are positioned 
symmetrically with respect to e g = 0. Furthermore the 
discrete energy spectrum shows the shell structure or the 
pair of tt and tt* levels. When the overlap between tt 
and tt* level is weak [panel (a) in Fig. 0] within one 
shell or one pair, the conductance peaks as well as the 
MR dips are distinguished between tt and tt* levels. The 
negative MR dips are clearly observed for each tt and 
tt* levels in the weakly overlapping case. With increased 
linewidth [panel (b) in Fig. @], the overlap between tt 
and tt* levels becomes significant and the MR structure 
is correspondingly modified. In our works, the Coulomb 
interaction between electrons are not taken into account. 
The Coulomb interaction, first of all, increases the level 
spacing between the conductance peaks so that the oscil- 
lating and negative MR structure will be persistent. 



Spin polarized transport through multi 
resonant levels 



Recently the research interest is growing in molecu- 
lar spintronic o 41 i 42 or spin polarized transport through 
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pends sensitively on the symmetry of the coupling matrix 
elements. Using the simple phenomenological model, we 
will show that the fine structure in magnetoresistance 
also depends on the symmetry in coupling matrix. 

The spin polarized current is computed using the same 
Eq. (jf 8|) . but the Green's function is now given in an en- 
ergy diagonal basis. Note that the Eq. (J20J) for the CNT 
Green's function is written down in the site-diagonal ba- 
sis of the tight-binding Hamiltonian. The two ferromag- 
netic electrodes are again modeled by the majority and 
minority spin bands, and the nanoscale system (NS) and 
the coupling between FM electrodes and NS are described 
by Hns and Hi , respectively. 



Ep/A 



Hns = EE E i4« d i<x, (28) 

i a— ± 

Hi = Yl T, V P^ k Hk a d ^+H.c (29) 
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FIG. 5: Effect of interfering energy levels on MR. A is the 
energy level spacing. Tl = ^Nl\Vl\' 2 and Tr are measured 
in unit of A. When the couplings to two ferromagnetic elec- 
trodes are symmetric and strong, neighboring energy levels 
overlap each other and the MR can become negative due to 
the interference between neighboring energy levels. Rather 
strong overlap of conductance peaks is needed to get the in- 
verse MR in the DBS. In the TSS, MR always remains to be 
positive in the symmetric coupling case. 



molecules. Our study in the CNT can be equally applied 
to the spin polarized transport through other nanoscale 
systems. Examples include the atomic wire, molecules, 
etc, which can accommodate the discrete energy spec- 
trum due to their finite size. The structure of the cou- 
pling matrix elements or the symmetry 4 ^ of wave func- 
tions in nanoscale systems plays an important role in 
determining the details of the electron transport through 
nanoscale systems. As an example, the existence of the 
transmission zeros (complete destructive interference) de- 



Ei = ei + Ep is the discrete energy spectrum of NS, which 
can be shifted by the gate voltage (Ep is proportional to 
the gate voltage). V P i is the coupling constant between 
the FM electrodes and the z-th energy level in NS. The 
retarded (r), advanced (a) Green's function of NS is 



D r f(e) = [e-E±iT a \ 



(30) 



where E^ — EiSij 



and T aij = ■nN La V iL V L j 



In our numerical simulation, we choose 



i x A and V iL = SiV L , V lR = Vr (real V L and 
Si = ±1 is the relative sign of the i-th energy level's 



nN Ra y iR Vnj 
e< = 
V R ) 

left and right coupling constants. The spin dependent 
linewidth is parameterized as T p ± = T p (l ± j3 v ), where 
(3 P is the spin polarization of the electrode p = L,R. 

The generic features of MR and conductance, which 
we found for the CNT, are also observed in the spin po- 
larized transport through multi resonant levels. Details 
won't be repeated here, but instead the effect of inter- 
fering conductance peaks on the MR will be discussed. 
In a resonant tunneling, there are typically two types 
of nanostructure o 44 ' 45 : double barrier structure (DBS) 
and t-stub structure (TSS). Though both DBS and TSS 
provide resonant energy levels for transport, the trans- 
mission coefficients are quite different due to the inter- 
ference between neighboring energy levels. TSS (sj = +1 
for all i) is featured with transmission zeros (complete de- 
structive interference) in every conductance valleys, while 
DBS [sj = (— l) 1 ] has no transmission zeros. As displayed 
in Fig. [5j the inverse or negative MR is possible in the 
DBS when the coupling to two FM electrodes are sym- 
metric and sufficiently strong. On the other hand the 
MR always remains to be positive or normal in the TSS 
when the couplings are symmetric, irrespective of their 
strength. In the case of asymmetric couplings, the MR in 
both DBS and TSS can be negative near the conductance 
peaks. 

The DBS is realized strictly in one-dimensional struc- 
ture like the carbon nanotubes, so that the inverse or 
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negative MR is possible in the strong symmetric cou- 
pling case. In the one-dimensional structure, the wave 
functions may have even or odd parity under the space 
inversion. This leads to alternating relative signs [si — 
(— l)*]i^ in the coupling matrix elements and the inverse 
MR in the strong symmetric coupling case. In the CNT, 
both 7r and n* states are one-dimensional so that the neg- 
ative MR is possible as shown in Fig. [3J The nanoscale 
systems can have the i-stub-like symmetry in their wave 
functions when they are extended in more than one space 
dimension. In summary, we found another case of the in- 
verse MR in the resonant transport when the couplings 
to FM electrodes are symmetric and rather strong. The 
interference between neighboring energy levels leads to 
the inverse MR 



IV. SUMMARY AND CONCLUSION 

Using the tight-binding Hamiltonian for tt electrons in 
a finite armchair CNT and the symmetry-adapted cou- 
pling constant between the FM electrode and CNT, we 
studied the conductance and MR in the spin valves with 
CNT. To characterize the CNT spin valves, we probed 
the model parameter space of the CNT Fermi energy, 
and the coupling strength between CNT and the two FM 
electrodes. 

In the case of asymmetric couplings between CNT and 
two FM electrodes, the MR has the broad positive peak 
at the conductance valleys, and the MR is dipped near 
the conductance peaks. Though the conductance shape 
is symmetric with respect to its peak position, the MR 
shape is asymmetric with respect to its dip position. 
The MR asymmetry is more enhanced with increasing 
coupling strength. The MR dip becomes negative with 
rather strong asymmetric couplings. The observed MR 
oscillation as well as the negative value o 14 ' 15 may well be 
explained by this model parameter regime. 

In the case of symmetric couplings between CNT and 
two FM electrodes, the MR is more sensitive to the cou- 
pling strength. In the weak coupling case, the MR is 
broadly peaked at the conductance valleys, and is posi- 
tive and suppressed with the dip-peak-dip structure near 
the conductance peaks. When the coupling strength is 
increased, the MR shape becomes highly asymmetric and 
the MR is negative near the conductance peaks. In 
Ref. ll_3j, the negative MR was observed for the highly 
transmissive contact between Co or Ni electrode and 
a CNT. In some of the spin valves, the CNT was to- 
tally submerged into the ferromagnetic electrodes. The 
highly transmissive contact means more or less symmet- 
ric couplings to the source and drain electrodes, and cor- 
responds to the case of a very strong coupling between a 
CNT and the FM electrode in our study. Our negative 
MR in the case of strong symmetric couplings might be 
able to explain the observed experimental results. 
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APPENDIX A: RELATIVE PHASE AT THE 
EDGE LAYERS: TIGHT-BINDING 
HAMILTONIAN FOR CNTS 

For completeness we include the tight-binding 
analysis 46 of a CNT electronic structure in this Ap- 
pendix, which is relevant to our study. The electronic 
structure of a graphene sheet can be successfully de- 
scribed by a single n orbital tight-binding model. As- 
suming the nearest neighbor hopping, the Schrodinger 
equation can be written as 

Ea = -t c j , ( A1 ) 

J6n.n. 

where t is the hopping integral which has the value of 
2.66 eV and Cj is the amplitude of tt orbital at the z-th 
site. The summation over j is restricted to the nearest 
neighbor sites of the site i. There are two atomic sites, 
C\ and C2, in a graphene unit cell (Fig. and the wave 
functions of neighboring unit cells are related to each 
other by the Bloch theorem. The Schrodinger equation 
can be written in a 2 x 2 matrix form using the Bloch 
theorem as 

<«) = -iV?)G)- (A2 » 

where k is the Bloch wave vector and h is given by the 
expression 

h(k x ,k y ) = 1 + e ils2 + e £ - Sl 

= l + 2cxpH^]cos^, (A3) 

with a — lai.,2 1 is the magnitude of lattice vectors which 
has the value of a = \/3 x 1.42 A ~ 2.46A. Referring 
to Fig. [HI note that ai, 2 = (±i + V3y)a/2. The energy 
bands E±(k) and the eigen states \E±(k)) are 

E±(k) = ±t\h(k)\, (A4a) 

\E±(k) > = -±= . (A4b) 

The phase 8 is defined by the relation e lS = h/\h\ and is 
given by the expression 

eiS = l + 2expH^]cosY (A5) 

^1 + 4008 cos Y + 4cos 2 %s 
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FIG. 6: Two atoms, Ci and C2, in a graphene unit cell and 
basis vectors, ai and a2. 



Now let us consider the energy bands in carbon nan- 
otube, which is formed by rolling up a graphene sheet. 
The circumference vector C'n,m — Na\ + MS2 uniquely 
dehnes the carbon nanotube. Periodic boundary condi- 
tion around the CNT circumference leads to the quan- 
tization of the Bloch wave vector along that direction 
or the distinct band indices. For the armchair type 
(n, n) tubes, C n .n — nv3a y and the quantization rule is 
k y = 2rmr /ny/3a. The (n,n) tubes repeat the unit cell 
(consisting of four atoms) n times around the circum- 
ference, so that there should be 4n energy bands in the 
(n, n) tubes. The band structure of (n, n) nanotube is 



E 



= ±^/l + 4cos^cos^ 



ka 
~2 



4cos 2 -^(A6) 



where m = 0, 1, ■ • • , 2n — 1. The relative phase difference 
between the wave functions of two atoms in a graphene 
unit cell can be written as 



l + 2exp[-i^]cos4f 



1 + 4cos^cos^ 



4 cos 2 !f 



(A7) 



For the case of m — and m = n, the band structure as 
well as the phase difference are simplified. 



E ±(k) = ±t 1 + 2 cos 



ka \ Ci 

, ka\ Ci 

B n± (k) = ±t ( l-2cos— J , -i 



T l,(A8a) 
T l.(A8b) 



Since we are interested in the spin polarized transport 
close to the Fermi level, we confine our interest to the 
m = n bands. 

Using the Bloch theorem, we can find the relative phase 
difference between two equivalent unit cells along the cir- 
cumference in (n, n) nanotubes, which is given by 



IV/) = exp[ik- (ai +a 2 )]\ip) 



. 2nm., , . , , „. 
exp[i M-(A9) 



Especially for the m = 0, n bands, electron's wave func- 
tions have a uniform phase from one unit cell to another 
around the CNT circumference. Combining all the rele- 
vant information, we summarize the m — 0, n bands and 
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(a) Schematic band structure 
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(b) Phase change around the circumference 

FIG. 7: (Color online) (a) Schematic display of 7r and 7r* 
bands. The lower ellipse represents the 7r band, while the 
upper ellipse does the 7r* band, (b) Relative phase between 
two atoms in a graphene unit cell. 



the relative phase in Fig. [7] The electronic states of the 
7r band (C1/C2 = 1) have uniform phase from one site to 
the next one in the edge layer. On the other hand, the 
7r* band has an alternating phase, ±1, from one site to 
the next one along the circumference direction. This dif- 
ference in phase will make the coupling strength between 
a CNT and the electrode depend on the symmetry states 
(tt or 7T*) in a CNT. 
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APPENDIX B: DETAILS ABOUT THE 
COUPLING BETWEEN CNT AND FM 
ELECTRODES 

In this Appendix let us study the coupling matrix be- 
tween the CNT and the electrode using the wave func- 
tions for free particles. We start with the wave functions 
for free particle in two-dimensional (2d) space. The ob- 
vious choice is the plane wave 



Here 



/V27T- For a free particle in 2d space, 
the wave functions for the 7r and 7r* states are given by 
the following 

{x\irkm) = \/2J nm (kr) cos (nm<^>), (BIO) 
(x\ir*km) = V2 Jnm(kr) svn{nm4>). (Bll) 

In general, for a given n, the wave functions can be writ- 
ten down as 



ipAx) = (x\k) 



2 7,2 



.4 



2m 



(Bl) 



Here A is the area of the two-dimensional system. The 
other one is the wave function with a cylindrical symme- 
try. 



^ km (x) = J m (kr)e m t 



2 1,2 



2m 



(B2) 



Here J m is the Bessel function. These two representa- 
tions of wave function for free particle in 2d space satisfy 
the orthonormality relation 



d 2 x ipl{x)ipp(x) = 



(B3) 



2tt 



d x ipl m {x)i>k>m>{x) = 5 mjm i— 8(k - k'). (B4) 



The following normalization is used for the cylindrical 
wave functions 

S e , e > = -^S(e-e') = ^(fc-fc') = 4. fc '.(B5) 
N 2 (e) k 

Here ^(e) is the density of states in 2d space and e = 
h 2 k 2 /2m and e' = h 2 k' 2 /2m. 

When the electrode is in contact with the armchair- 
type (n, n) CNT, electrons in the electrode should belong 
to the same symmetry states as in the CNT if they have 
any chance to hop into the CNT. Since only the n and 
it* states are responsible for the electron transport near 
the Fermi level under the neutral charge condition, we 
consider only the irreducible wave functions belonging to 
the 7r and it* states. For the (n, n) CNT, 



1 



7T 777) 



COs(77777(/)), 



: Sin(777710). 



(B6) 
(B7) 



(x\Trkm) = R m (kr) cos (nmq 
(x\n*km) — R m (kr) sm(nm4 



(B12) 
(B13) 



That is, the angular dependence is determined by the ax- 
ial symmetry and does not change for any realistic wave 
functions. However, the radial part, reflecting the com- 
plexity of the system, will be more complicated than the 
Bessel function. 

Let us consider the coupling between the (n, n) CNT 
and the ferromagnetic electrode. Using the irreducible 
representations, the coupling Hamiltonian 38 can be ex- 
panded as 

Hi = Y, c Uk m (^ km W\j)a J +H.c. 

qkm 



+ Y, c l* q k m (**<lkm\V\j)a J + H, 

qkin 



(B14) 



The wave vector K — qz + k is defined in the electrode, q 
is the wave number along the CNT axis and k is normal 
to the CNT axis. The omitted part (neither tt nor it* 
states) is not relevant to transport near the Fermi level, 
and will be neglected from our study. In the tight-binding 
Hamiltonian approach, we adopt the localized Wanier- 
type wave functions for the p z electrons at the carbon 
site. The coupling matrix can be written down as 

(■jrqkm\V\j) = V„ m (q, k) cos(nm^), (B15) 
(ir*qkm\V\j) — V n * m (q, k) sm(nm(pj). (B16) 

Here <pj represents the angular position of j-ih carbon 
atom along the circumference at the edge layer of a (77, 71) 
CNT. In our study V^ m {<l,k) and 14-* m (g, A;) are treated 
as phenomenological parameters. When the angle <pj is 
substituted into the above coupling matrices, we get the 
desired simple forms 



Here <f> is the azimuthal angle around the CNT axis. From 
these irreducible wave functions, we can deduce the fol- 
lowing identities (777 = — m) 



|7T777 > = ^p(|?777l) + |n.T7l)), 

v2 



7r m > 



V2 



(— \nm) + \nm)). 



(B8) 
(B9) 



777 7T 



(Trqkm\V\j) = V vm (q, k) cos — , 



(B17) 



(tt* qkm\V\j) 



(-l) J 'K. m (g,fc)sin— . (B18) 



The factors from angular functions are included in the 
coupling strength. As expected, the relative phase along 
the circumference of one layer correctly matches the it 
and 7r* states in the CNT. 
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TABLE I: Character table for Cs v and decomposition of re- 
ducible representations, a = 2cos27r/5 = (—1 + \/5)/2 and 
b = 2cos47r/5 = (—1 — Vo)/2. 7^: representation of ir or- 
bitals. 7 m : representation of angular momentum states m 
and — m. 





E 


2C 5 


2Cl 


5a v 


A 1 


1 


1 


1 


1 


A 2 


1 


1 


1 


-1 


Ex 


2 


a 


b 





E2 


2 


b 


a 







20 











70 


1 


1 


1 


1 


7™^o 


2 


2C0S^2L 


2cos^ 

O 






APPENDIX C: SYMMETRY IN COUPLING 
MATRIX ELEMENTS BETWEEN A CARBON 
NANOTUBE AND FERROMAGNETIC 
ELECTRODES 

Let us study the symmetry in the coupling matrix be- 
tween the carbon nanotube and the ferromagnetic elec- 
trodes. Group theory^ is quite a powerful tool for the 
discussion about the symmetry in the coupling matrix 
elements^ Since the armchair CNT has the discrete ro- 
tational symmetry about the carbon nanotube axis, the 
energy bands or levels can be classified according to the 
irreducible representations of the symmetry group of a 
CNT. Correspondingly the conduction electron states in 
the FM electrodes should have the same symmetry as the 
CNT states in order for electrons to be able to traverse 
the CNT. 

As a concrete example, we consider the (5, 5) armchair- 
type CNT with a finite number of layers. Our discus- 
sion can be easily extended to the general (n, n) arm- 
chair CNT. Depending on the even or odd layers, a finite 
CNT may have an additional symmetry along the nan- 
otube axis. Since the symmetry in the coupling matrix 
is mainly determined by the rotational symmetry about 
the CNT axis, we are going to confine our discussion to 
this rotational symmetry. 

The finite (n, n) CNT belongs to C nv symmetry 
group4£ For the (5,5) CNT, there are four distinct 
classes: E (identity), 2C5,2C| (five-fold rotations) and 
5a v (reflections). Correspondingly there are four irre- 
ducible representations for this group. The character ta- 
ble is summarized in Table [J and can be used to find the 
irreducible representations for any reducible representa- 
tion. 

First we consider the electronic states in an infinite 
CNT. For the electronic structure, we adopt a single tt 
orbital approximation. In this case the p z orbital is lo- 
cated at each carbon site. Acting the symmetry opera- 
tions on the CNT, we can find the characters of the 7r 
orbitals. There are 20 carbon atoms for (5, 5) CNT in an 
extended unit cell around the circumference. Using the 
decomposition formula, we can readily find the number 



of irreducible representations. 

7 T = 2A-L + 2A 2 + 4Ei + 4B 2 . (CI) 

There are four nondegenerate bands and eight bands with 
double degeneracy. When we take into account the phase 
variation of the tt and tt* states along the circumference 
of the CNT, we can deduce that tt states belong to the Ai 
irred. reps. On the other hand, the tt* states belong to 
the A2 irred. rep. Note that 7r and tt* states are invariant 
under discrete C5 rotations. The tt states are invariant 
under a v , but the tt* states change the sign under a v . 

Now let us consider the angular momentum states 
about the CNT axis. The angular momentum operator 
L z is defined about the nanotube axis and its eigen states 
are well known to be ipm{(f>) = e w ™* > , with L z ip m ((f)) = 
mip m (4>). Since ip m is the basis function for the contin- 
uous rotational symmetry, it cannot be the eigenstate 
for the discrete rotational symmetry, e.g., C§ v for (5,5) 
CNT. Acting the symmetry operations on these basis 
functions, we can find their character tables. The effect of 
the identity E, and the discrete rotations C5, Cf , Cf , C| 
on the angular momentum wave functions is obvious. 

E^ m {4>) = VwW, (C2a) 

K a ^ m {<\>) = VmW-a) = e~ ima ^ m (<t>) . (C2b) 

The second operation TZ a is to rotate the wave function 
by the angle a. The characters can be easily identified. 
For the operation of a v , a reflection, we have to find out 
its effect on the angular momentum operator L z . The ef- 
fect of <j zx (zx plane) and a zy (zy plane)on L z is obvious. 

o zx L z a^ x = —L z , (C3a) 
G zy L z a~y = -L z . (C3b) 

Note that x — > x,y — > —y,z — > z and p x — > p x ,P y — > 
—PyiPz —> Pz under a zx . Let us see if the above relation 
is true for an arbitrary reflection a v where the z axis lies 
on the reflection plane. 

o v L z g~ x = -L z l (C4) 

We can prove this relation by the direct construction. 
Assume that the angle between the plane a v and the 
plane a zx is a. Under a v , the coordinate (x, y) transforms 
into a new one (a/, y'). They are related by the equations 

x' — £ cos 2a + y sin2cv, (C5a) 
y = x sin 2a — y cos 2a. (C5b) 

The linear momentum operators transform in the same 
way. Then wg CEtn readily show that L z — %py — UPx 
transforms as 

<J. v L z o- x = -L z . (C6) 

Now let's see the effect of a v on the wave function ip m {(/>)- 
Intuitively o v changes the direction of a rotation: right- 
hand rotation into left-hand rotation. Hence we can de- 
duce that 

cr v ip m ((j)) = ip m (-(f>) = ip-m(<f>). (C7) 
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We can prove the above relation in a formal way using 
the transformation rule for L z . Consider a v L z ip m (<p) = 
fn<^v'4 ) m{4 1 )- On the other hand, a v L z ip m ((f)) = 
o v L z o~ x o v ^ rn {4>) = -L z a v Tp m (<fi). Since L z a v ip m (4>) = 
—mo v ip m {4>), we can deduce that 

0- v ^m{4>) = c i>-m(<f>)- (C8) 

Here c is the C-number to be determined. Double oper- 
ation of a v is an identity so that c? = 1. We may choose 
c = 1. Since a v couples two states, </> m and ip-m, we have 
to consider two states on the same footing. Let j m denote 
the (reducible) representation in the functional space of 
(ipmjip-m)- We can now build up the character tables 
for 7 m considering the symmetry operations on the basis 
functions. For example, 

'.(&)-(:s)(&)- <»> 

Hence the character Xm(<^v) = (the trace of a transfor- 



mation matrix). The other characters can be obtained 
using the same process and summarized in the Table HI 
Using the decomposition rule, we can identify the irre- 
ducible representations. 



To = Ax, (ClOa) 

75m = A x +A 2 , m^O, (ClOb) 

75m±l = Ex, (ClOc) 

75m±2 = E 2 . (ClOd) 



From the above group theoretical analysis, we conclude 
that the it states in a finite (5,5) CNT are coupled to the 
angular momentum states cos 5m0, while the n* states 
are coupled to the states sin 5m0. Other angular momen- 
tum states are coupled exclusively to the states belonging 
to the Ex and E 2 irreducible representations, but never 
coupled to the it and ir* states. 
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